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Abstract — Most cellular radio systems provide for the use of 
transmitter power control to reduce cochannel interference for a 
given channel allocation. Efficient interference management alms 
at achieving acceptable carrier-to-interference (C/I) ratios in all 
active communication links in the system. In this paper, such 
schemes for the control of cochannel interference are investigated. 
The effect of adjacent channel interference is neglected. As a 
performance measure, the interference (outage) probability is 
used, Le^ the probability that a randomly chosen link is subject 
to excessive interference. In order to derive upper performance 
bounds for transmitter power control schemes, algorithms are 
suggested that are optimum in the. sense that the interference 
probability is minimized. Numerical results indicate that these 
upper bounds exceed the performance of conventional systems 
by an order of magnitude regarding interference suppression and 
by a factor of 3 to 4 regarding the system capacity. The structure 
of the optimum algorithm shows that efficient power control and 
dynamic channel assignment algorithms are closely related. 



L Introduction 

EFFICIENT channel reuse is of paramount importance 
in the design of high capacity cellular radio systems. 
Cochannel interference caused by frequency reuse is the 
single most restraining factor on the system capacity. Systems 
may be designed either to tolerate low signal-to-interference 
ratios (Le., efficient modulation and coding schemes), or 
to reduce the interference with efficient cell planning. In 
addition, dynamic channel allocation may be used to reduce the 
interference even more by adapting the cell plan to changing 
traffic and unforeseen propagation conditions. 

In addition, the use of transmitter power control has been 
proposed to control cochannel interference. The main idea is 
to adjust the power of each transmitter for a given channel 
allocation, such that the interference levels at the receiver 
locations are minimized. Maintaining adequate transmission 
quality on the actual communications links is an obvious 
constraint The measure of quality usually employed in cellular 
systems design is me caaiet^^ 

benefits of power control are, however, not as obvious as they 
may seem at first glance. It should be borne in mind that by 
reducing the transmitter power in a certain link, that link will 
also be more vulnerable to interference. That the design of 
power control algorithms is not trivial is illustrated by the 
wide range of algorithms and results found in the literature. 
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Several aspects of power control for cochannel interference 
reduction have been investigated in the literature. Early an- 
alytical work by Bock and Ebstein [1] has developed into a 
well-known technique in government spectrum management 
[17]. The problem is here approached by means of linear 
programming. In [2], Aein investigates cochannel interference 
management in satellite systems. He introduces the concept 
of C/I balancing, which yields a "fair" distribution of the 
interference in the sense that all users experience the same 
C/I level. The problem is identified as an eigenvalue problem 
for positive matrices. In [3]-[5], Nettleton and Alavi apply 
and extend these results to spread spectrum cellular radio 
systems. In these systems, the adjacent channel ("adjacent 
code") interference also has to be taken into account. In fact, 
this type of interference is shown to be dominant in spread 
spectrum systems [18]. Nettleton and Alavi show that C/I 
balancing substantially improves the capacity of such systems. 

Other power control algorithms implemented this far basi- 
cally keep the received power at a constant level This scheme 
is very efficient in controlling adjacent channel (intracell) 
interference, but has been shown to have limited effect on 
cochannel interference [9]. Most work done in this area 
involves simulation studies with quite realistic but complex 
system models {6]-[10]. Unfortunately, the large variety of 
model assumptions makes these results hard to compare. 

The aim of this paper is to further investigate the per- 
formance of transmitter power control algorithms. The basic 
models and ideas in [2]-[5] form the point of departure for 
the following discussion. Hie approach in this paper will be 
to find performance bounds and conditions of stability for all 
types of transmitter power control algorithms. To derive these 
bounds, a (hypothetical) power control scheme is proposed 
that is optimum in the sense that it ininimizes the interference 
(or outage) probability. The latter quantity is the probability 
that some r andomly chos en mobile (or base station) has aJ3ta 
low" C/I. It will be seen that the C/I balancing technique 
plays a key role in this minimization process. Finally, some 
implementational aspects regarding distributed C/I control 
algorithms will be discussed. Hist, however, we will define 
our system model and derive some useful results from the 
theory of positive matrices. 

n. System Model 

Throughout this paper we will study a large, but finite, 
cellular radio system consisting of N cells. For the system 
design we have at our disposal M independent channel pairs, 
each consisting of independent up- and downlink channels. 
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This means that the crosstalk between channels (adjacent 
channel interference) is neglected. The actual implementation 
of the channel is not of importance to our model. Frequency 
slots (FDMA), times slots (TDMA), or even spreading codes 
(CDMA) may be used. In CDMA systems, however, the 
assumption of having no crosstalk between channels may 
not be very realistic. The channel pairs will be allocated to 
mobiles and base stations according to some arbitrary channel 
assignment strategy. Each channel pair is assumed to be used 
only once in each cell, which means that there will be no 
intracell interference. The set of cells using a certain channel 
pair m at some given instant, will be called the cochannel set 
of m. The size (cardinality) of this set is denoted Q(m). 

Let us now consider a call in progress in some cell of the 
cochannel set of m. Both the base station and the mobile 
station will experience interference from the other Q(m) — 1 
cochannel cells of the set We will in the sequel assume the 
transmission quality to be dependent only on the C/I ratio, 
I\ The total interference power is modeled as the sum of 
the powers of all active interferers [13], [16]. We will further 
consider high capacity, interference limited systems only. Hie 
consequence of this is that other noise sources, like termal 
noise, may be neglected. Using these assumptions we may 
express Y as 

r = Prs%i (1) 

Q(m)-1 V 9 

where P n ,i is the received power from "the "desired" trans- 
mitter and Ij is the received power from interferer j. 

Now, let us denote the link gain on the path between the 
mobile in cell i and the base station in cell j at some given 
moment by Gij (Fig. 1). Note that, due to our assumptions 
above, there is only one mobile in cell i using this particular 
channel. Using this notation, we may derive the C/I at mobile 
i 

P rx , GaPi 



CcUi 



HI, 



3& 



(2) 



where Pi is the transmitter power used by the base station in 
cell i. Note that Gu is the path gain of the "desired" signal path 
in cell i. To simplify further derivations, we introduce some 
matrix notation similar to that in [2]-[5]. First, we define the 
link gain matrix G = {G^}. It should been borne in mind 
that, in a mobile system, the link gains will constantly change. 
G may, therefore, be described as a matrix valued random 
process. At some given instant of time, however, G will be 
a Q x Q square matrix of random variables. Further dividing 
the right-hand side (RHS) of (2) by. Gu yields 



(3) 



J5« 



i=i 

where we have introduced the stochastic variables Zij defined 
as 




Fig. 1. Link geometry and link gain model. 

Zij is the link gain from the base station in cell j to the mobile 
in cell i, normalized to the link gain in the desired path in cell 
i (i.e., from base station i to mobile i). Again using matrix 
notation we introduce the normalized downlink gain matrix 
Z — {Zij}. We may in the same manner derive the C/I at the 
base station by letting P denote the mobile transmitter powers, 
and replacing Z with the uplink gain matrix W — {Wij} 
defined as 



Wi 



Gu „ Gu 



J J3 



Zij = 



(4) 



Note that in general Gu is not equal to Gjj. Therefore, W 
will, in general, not be identical to the transpose of Z\ The 
properties of a given uplink matrix are thus not identical to 
the properties of the corresponding downlink matrix. However, 
it is reasonable to assume that they exhibit the same statistical 
properties. 

Ti as defined by (3) will be a stochastic variable. Now, let 
T denote the C/I at some randomly chosen mobile. We define 
the distribution of T as 

l Q 1 Q 

F< 7 ) = Pr{r < 7} = 5 E Pr ( r i * 7} = ^ X>(7) 

(5) 

where Fj(i) is the distribution of Tj. In general, the distribu- 
tions FjQy) are not identical for all j's. In all finite networks 
there will be boundary effects, slightly favoring boundary cells. 
Let us now assume that the transmission system requires a 
minimum (threshold) C/I ratio 7c This C/I level is called 
the system protection ratio. As our performance measure, we 
will throughout the paper use the interference probability 

F( 7 o) = Pr{r<7o}- (6) 

This is the probability that some randomly chosen mobile has 
a C/I below the system protection ratio [13], [16]. 

Now, let us define a power control algorithm (PCA) as an 
algorithm using some measurement information to determine 
the transmitter power vector P = {Pj}. The performance 
of the algorithm is measured by the resulting interference 
4 probability, Ffro). We will define a global PCA, as 
an algorithm that in every moment has access to the entire 
link gain matrix Z and may instantaneously control the entire 
vector P. can be described by a vector valued function 

P = V G (Z) (7) 
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In a local PCA, however, each base station (or mobile) controls 
its own transmitter power based on only limited knowledge 
about Z. It is clear that the class of local PCA's is included in 
the class of global PCA's. Therefore, the lowest interference 
probability that may be achieved by a local PCA cannot be 
better than the interference probability achieved by the best 
global PCA. The global algorithm may, therefore, be used 
to derive an upper bound on the performance of any local 
algorithm. From a practical point of view, however, local 
algorithms are preferred since they are simpler and require 
considerably less information exchange. 

III. Achievable Cochannel Interference 

Let us consider a system with global power control and 
some given instantaneous link gain matrix Z, 

Definition: The C/I-levei 7 is achievable if there exists a 
power vector P > 0 such that r» > 7 for all cells (mobiles) i. 

The notation P > 0 is used to denote that all elements of 
P are nonnegative. The following result is extremely useful 
to determine the performance of a global PCA. Given a 
nondegenerate stochastic link gain matrix Z we have the 
following proposition. 

Proposition 1: With probability one, there exists a unique 
maximum achievable C/I level 

7* = max{7|3 P>0:Fi> 7,V £}. 

The maximum is given by 




where A* is the largest real eigenvalue of the matrix Z. The 
power vector P* achieving this maximum is the eigenvector 
corresponding to A*. 

Proof: The matrix Z is by definition a nonnegative 
matrix (all elements are nonnegative). Since all elements of 
Z are stochastic variables, the matrix will (for nondegen- 
erate cases) have full rank with probability one. From the 
theory of nonnegative matrices we collect the following results 
originally due to Perron, Frobenius, and Wielandt [11]. 

i) Z has exactly one real positive eigenvalue A* for 
which the corresponding eigenvector is positive (i.e., 
all components have the same sign), 
ii) The minimum real A such that the inequality 

AF > ZP (8) 

which has solutions for P > 0 is A = A*. 
Now, rewriting (3), our proposition states that 

i=i 

or 

^Pi>J2PjZ^ Vi. 
7 j=i 



Rewriting this in matrix form yields 

±±^P > ZP. 
7 

Combining results i) and ii) yields the required result. The 
simple bounds on the eigenvalues of Z derived from row sums 
of Z given in [2], [11], ensure that A* > 1 and thus 7* > 0. 

□ 

Proposition 1 tells us that the maximum achievable C/I 
level is in fact reached when all inequalities (9) are satisfied 
with equality. This means that a C/I balanced system as 
proposed in [2]-[5] reaches the largest achievable C/I level, 

IV. Optimum Power Control^ 

Now, let us assume that we use a global PCA, i.e., a vector 
valued function defined in (7). Let the achieved C/I vector 
given in (3) be denoted IT = {ri} where 

r = r(p, z) = it(#g(s), z). (io) 

Now using (5) and conditioning on the link gain matrix, Z 9 
we have 

1 Q 

F(7o|^; *g) = 75 £^(7012; **)- (11) 

Now, since &g is a deterministic function and Z is given, the 
distribution Fj becomes a simple step-function, i.e. 

F (^\Z* - / ll r ' ( * G ' Z) < To 

F wo\ z > * G > - \ 0, T^g, Z) > 70. C } 

We may thus express (11) as 

F(7o|2,*g)= ±B{Z,* g ) (13) 

where the integer function B is defined as 

B{Z, ^ G {Z)) = # elements mT(^ G (Z)Z) < 7o. (14) 
Taking the expectation over Z we obtain 

F(rro;*G) = ^E[B{Z,* G )] 

= = (15) 

To minimize -F(7o) we seek to minimize the function B{Z^g) 
for every outcome Z. This minimum number we denote by 
B*{Z). B*(Z) is the minimum number of inequalities (9) that 
do not hold, ie., the number of links in which the protection 
ratio 70 may not be achieved. Since B* is an integer function, 
the minimum may be reached by several transmitter power 
vectors. In particular, we have the following proposition. 

Proposition 2: At least one optimum power vector P has 
the form 

Pi = 0, i € sn 
where SH = {% : Ti < 70}, |«| = B*. 
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Proof: Assume wc have an optimum power vector P* 
not fulfilling the above requirement, Lc, P) £ 0 for all j. Let 
r* be the corresponding C/I vector. We form a vector P such 
that, for some k € 9t : 

Pk = 0. 

Computing the C/I vector T for this power vector we obtain 

r i= — * = _ 2 

£ Pi 2* - ^ E *7 ' ^ - p * z " 

r fc = o 

since an quantities are nonnegative. P achieves a C/I that is 
not less than 70 whenever this is achieved by P*. Therefore, P 
is also an optimal power vector. Lowering the C/I for the link 
in cell k will be of no consequence since the link is considered 
to be unusable anyway. Repeating the argument stepwise for 
all indexes k € 9* proves the proposition. D 
Due to Proposition 2, we may limit our search for optimum 
power vectors to those vectors with positive components for 
those cells where the required C/I is achieved. The other 
components are zero, corresponding to no transmission at all. 
By setting some P k = 0, we effectively remove cell k from 
the cochannel set. Now, let us form a submatrix of Z, denoted 
Z* where all rows and columns k corresponding to zero 
components in the optimum power vector have been removed. 
Clearly, the maximum achievable C/I-level for Z* must be 

not less than 70- . . . 

This observation leads us to a procedure to minimize the 
outage probability. An optimum global PCA will, by removing 
as few cells as possible, find the largest submatrix Z* for 
which 70 is achievable. In this way we find an optimum 
power vector consisting of the eigenvector of Z* and zero 
components corresponding to the removed cells. The "cell 
removal" procedure may appear as a strange way to achieve 
high system capacity. However, one should bear in mind that 
we consider those links* where the required C/I level cannot 
be reached to be useless. Using a positive transmitter power 
in such a link will only cause interference in other cells. Note 
that straightforward C/I balancing, without cell removal, may 
be disastrous since all links may drop below the C/I threshold. 

The question arises how one would implement such an 
optimum procedure. A "brute force" search would first check 
if 70 is achievable for the original matrix Z. If not, we would 
try to remove one cell, computing the eigenvalue of each 
reduced system until the C/I requirement was fulfilled. If not, 
we would try removing all combinations of 2 cells and so 
on. Clearly, such a procedure would find an optimum power 
vector in a finite number of steps since finally, by removing 
combinations of Q - 1 cells, no interference would remain. 
However, in the "worst case," the total number of eigenvalue 
computations we have to perform is exponentially increasing 
with Q- To overcome this practical problem we propose the 
following algorithm that in most cases closely approximates 
the behavior of the optimum algorithm. We use a sequential 



procedure, where each step involves only one eigenvalue 
computation. The idea is to remove one cell at the time until 
the required C/I level is achieved in the remaining cells. An 
obvious criterion for cell removal is of course to maximize the 
(7')* of the reduced system. This, however, requires Q-k 
eigenvalue computations at the fcth algorithm step (and roughly 
a total of Q 2 /2 computions in the worst case). Instead, we 
propose an even simpler criterion that requires substantially 
less computational effort. 

Stepwise Removal Algorithm (SRA) 

1) Determine 7* corresponding Z. If 7* > 7o use the 
eigenvector P*, else set Q' = Q. 

2) Remove the cell k for which the maximum of the row 
and column sums 

Q 

r fc =VZ*i, r% = 22Z jk 

is maximized (combined sum criterion) and form the (Q' - 
x \ x (Q' _ 1) matrix Z'. Determine 7* corresponding to Z . 
If 7* > 7o, use the eigenvector P*, else set Q' = Q' - 1 and 
repeat step 2). 

By this procedure we, one by one, remove cells stepwise 
until all the C/I's in all remaining cells are larger than 70- 
The row and column sums provide bounds on the dominant 
eigenvalue of the matrix Z [2], [11]. The cell removal cri- 
terion, therefore, seeks to maximize the lower bound for the 
7 * of the next matrix Z\ Other simple criteria have been 

evaluated in [15]. 

At the present we do not know if an algonthm exists to 
solve the original submatrix problem in polynomial time or 
if this problem is actually JVP-complete. Thus, we cannot 
determine if, and under what conditions, the problem in 
fact may be solved by a step-by-step procedure like the 
SRA. This constitutes an interesting open problem for further 
investigations. 



V. Numerical Results 
In order to derive some numerical results, more specific 
model assumptions have to be made. We will, for this purpose, 
consider a fixed and symmetric channel assignment strategy 
that divides the cells in K different channel groups. The cells 
using the same frequency are place symmetrically in a hexag- 
onal grid [12]. Each set will consist of N/K ceUs forming 
a sparse hexagonal pattern. Base stations use omnidirectional 
antennas and are located at the center of the cells. The locations 
of the mobiles are assumed to be uniformly distributed over the 
cell area. To simplify calculations we will express all distances 
normalized to the radius of the hexagonal cells. 

A simple propagation model that has found widespread use 
in the analysis of cellular systems is used [12]. We will assume 
that adequate modulation, coding, and equalization schemes 
are used in order to suppress the explicit dependence of the 
multipath fading process. Therefore, the performance of the 
system will only be dependent on the "local average" received 
signal power. The link gain G i5 is modeled as 



!.. = Ail 



x *3. 



(16) 
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Fig. Z Distribution of the maximum achievable C/I in all cells in the 
.system, 7*(k), when k cells are removed according to the combined sum 
criterion. The system investigated contains 16 cochannel cells -and uses a 
seven-cell cluster. Log-normal fading, <r = 6 dB, a = 4. 

where d+j is the distance between the mobile in cell i and the 
base station in cell j. The 1/cP factor models the large scale 
propagation loss whereas the attenuation factors Aij model 
the power variation due to shadowing. We assume all A*j 
to be indentically independent and log-normally distributed 
random variables with 0 dB expectation. The log-variance will 
be denoted by a. Parameter values of a in the range of 3-5 
and a in the range of 4-10 dB usually provide good models 
for urban propagation [12]. 

The numerical results are for a 16 cell "square" cochannel 
cell pattern with "full load," Le., all cochannel cells are 
assumed to be in use. The evaluation of the interference 
probabilities is made by means of Monte Carlo simulation. 
The standard deviation of the estimated probability at the 10% 
interference probability level is around ±1%. 

Due to the excessive computational complexity of the "brute 
force" optimum algorithm this paper is mainly limited to the 
stepwise removal algorithm. Unless otherwise noted, we use 
the propagation constant a = 4 in the following graphs. Lower 
values of a yields somewhat more pessimistic interference 
probability results. However, the relative performance of the 
different algorithms investigated hardly changes at all. 

Fig. 2 gives an insight into the operation of the SRA 
algorithm. It shows the achievable C/I levels, *y*, in the system 
where k cells have been "removed" by the SRA. As can 
be seen from the graph, the expected achievable C/I level 
increases. rapidly as a few cells are removed according to 
the removal criterion. Apparently very few cells cause the 
most interference problems. These cells are often, but not 
- always, "interior" cells of the cochannel cell pattern. Results 
further indicate that, for a low number of cell removals, the 
performance of the SRA algorithm is very close to optimum. 

Figs. 3. and 4 show the performance (interference proba- 
bility) of the SRA algorithm compared with the performance 
of two reference algorithms, fixed transmitter power (Lei, no 
power control at all) and constant received power. The cluster 
sizes K = 3, 7 are covered in these graphs. Results show that 
the SRA algorithm outperforms the reference algorithms by 
an order of magnitude. At the 10% outage probability level 
C/I gains in excess of 10 dB may be achieved. The fixed 
transmitter power and the constant received power algorithm 
achieve roughly the same performance for interesting outage 
probabilities. 



Fig. 3. Interference probability comparison. The system investigated con- 
tains 16 cochannel cells and uses a three-cell cluster. Log-normal fading, 
o = 6 dB, a = 4. 




Fig. 4. Interference probability comparison. The system investigated con- 
tains 16 cochannel cells and uses a seven-cell cluster. Log-normal fading, 
a = 6 dB, or = 4. 




Fig. 5. Interference probability for different cluster sizes. The system inves- 
tigated contains 16 cochannel cells and uses a 3, 7, 9, and 13 cell clusters. 
Log-normal fading, a — 6 dB, or =s 4. 



Fig. 5 finally, shows a comparison between-systems of 
different cluster sizes. It can be seen clearly that a system 
with cluster size 3 using the SRA algorithm outperforms a 
cluster size 13.system using the reference algorithms for most 
protection ratios. Capacity gains in excess of 4 times are, 
therefore, to be expected. 

VL Discussion 

In this paper we have studied performance bounds for 
power control algorithms. Our approach has been to study the 
general properties of some hypothetical global power control 
algorithms (PCA's). These algorithms assume knowledge of 
the attenuation in all transmission and interference paths. A 
PCA is described that is optimum in the sense that it minimizes 
the interference probability for a given threshold C/I (i.e., 
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protection ratio). This algorithm was, however, found to be 
somewhat impractical, since the computational complexity 
was exponentially increasing with the size of the cochannel 
set However, a sequential algorithm that approximates the 
behavior of the optimum algorithm in "linear" time, the SRA, 
is presented. The numerical results presented for this algorithm 
indicate that there are potentially large gains to be achieved by 
using C/I control Gains in excess of 10 dB in interference 
reduction may be achieved. This would correspond to capacity 
gains in the order of a factor 4 compared to system using fixed 
transmitter power or traditional power control. However most 
of the derivations are concerned with the downlink channel, 
results are equally valid for the uplink channel. It should 
be noted that the algorithms discussed in Section IV are not 
dependent on any of the explicit model assumptions used for 
deriving the numerical results. 

All these algorithms apparently involve quite drastic mea- 
sures to achieve a low interference probability. Usually, in 
one or more cells, base stations and/or mobiles are required to 
discontinue their transmissions. This is, however, a straight- 
forward consequence of the performance measure. Links that 
end up with a C/I level below the required protection ratio 
are considered useless. In digital communications with exten- 
sive error control coding this is a quite realistic assumption. 
Wasting transmitter power in these cells will only add to 
the interference in the system. In fact, this is the usual 
remedy to poor C/I levels in modern cellular systems. A 
mobile experiencing poor transmission quality is handed over 
to another channel and/or base station which improves the 
performance locally. The power control algorithms discussed 
above tell us which of the mobiles that should be handed 
over to maximize global system performance. The intimate 
relation between dynamic channel assignment and efficient 
power control is obvious. 

Attempting a practical implementation of the SRA algorithm 
would be difficult, since the path gain elements of the (time 
varying) link gain matrix Z, are usually not known. Estimating 
these gains would require a significant measurement effort. 
Even if this would be possible, the amount of data that 
would have to be communicated and managed by the central 
controller would be enormous in a reasonably sized network. 
The main benefit of these results is, therefore, that they provide 
an estimate of the optimum performance of all C/I control 
schemes. 

Further studies should, therefore, be devpted to studying 
algorithms that use only limited path gain information. Of 
special interest are local PCA's that may be operated only with 
path gain information that can be derived from measurements 
in a single cell. One interesting direction would be to inves- 
tigate distributed C/I balancing algorithms, i.e., algorithms 
that attempt to achieve a given C/I level [14]. The result in 
Proposition 1 would establish for which preset C/I levels such 
an algorithm would be stable, Le., that C/I balance could be 
reached. 
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